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Abstract 

We introduce a fundamental theory for the kinetics of systems of classical particles. 
The theory represents a unification of kinetic theory, Brownian motion and field theory. 
It is self-consistent and is the dynamic generalization of the functional theory of fluids 
in equilibrium. This gives one a powerful tool for investigating the existence of ergodic- 
nonergodic transitions near the liquid-glass transition. 
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I. INTRODUCTION 

We present here a theory^ for the dynamics of classical particles which solves the chronic 
problem of self-consistency. This theory unites the desirable elements of kinetic theory^, 
Brownian motion^, and modern field theory^. 

Kinetic Theory is one of our oldest^ theoretical disciplines. Despite its many successes 
it has never been constructed in a fully self-consistent form^. Thus one of the most famous 
approximations in all of science, the stosszahlansatz of Boltzmann^ and the treatment of 
the collision integral in the Boltzmann equation^ have not been investigated systematically. 
The theory introduced here provides the tools to remedy this situation. 

The stosszahlansatz, also referred to as the assumption of molecular chaos, is repre- 
sentative of decoupling approximations appearing in manj^ problems and characterized as 
uncontrolled by the approximation police. Of particular current interest is the validity of 
mode coupling theory (MCT)2 used in theories of the liquid glass transition™. We should 
be able to answer the question: Is the liquid-glass transition accompanied by an ergodic- 
nonergodic (ENE) transition?— The construction^^ of mode coupling models using tradi- 
tional kinetic theory is ad hoc and short of convincing. It has been completely ineffective 
in exploring corrections to conventional mode coupling theory. The field has moved away 
from kinetic theory treatments and turned instead to field theoretical models^^ where one 
has the promise of perturbative control. Thus we recently introduced the random diffusion 
model^ which can support an ergodic-nonergodic transition^ at one loop order. However 
going to two-loop order one finds that the system can not sustain the ENE solution. The 
problem with such models, compared to microscopic models organized in terms of a pair 
potential, is that the short- distance structure is not treated naturally. In the case of col- 
loidal systems, there has been some convergence on the Dean-Kawasaki (DK) model^ as the 
simplest field theoretic model that describes the kinetics of the colloidal systems operating 
under Smoluchowski^ dynamics. It has been difficult, for technical reasons, to establish 
whether the Dean-Kawasaki model supports an ENE transition even at one loop order. 

^ The title of this paper suggests that the theory is more widely applicable than to Smoluchowski dynamics. 
This is true. The basic development goes through for Newtonian dynamics, Fokker-Planck dynamics 
and trapped systems. However, presenting the theory in its most general form makes the development 
more complicated. Instead, I have chosen to present the theory in its simplest application, the case of 
Smoluchowski dynamics. The case of Newtonian dynamics will be available soon (S. Das and G. Mazenko). 
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It is demonstrated here how mode couphng theoryi^ naturally occurs in the application of 
our theory to colloidal systems governed by Smoluchowski dynamics. The approach, which 
allows for compatible approximations for higher order correlation functions^, is applicable 
to a large set of dynamical systems,^- reversible^ and dissipative, including Newtonian^!, 
Fokker-Planck=2-- and Smoluchowski^^ dynamics. The theory is organized in terms of a cou- 
pling to time and space dependent external fields. This allows for great flexibility in using 
functional methods in developing various types of perturbation theory. 

A key point is that the equilibrium equal-time fluid structure^^ has been understood 
from a self-consistent field theoretical point of view for a long-time. The work presented 
here is the natural extension to the dynamic regime of the beautiful diagrammatic/functional 
development for the static properties. The theory in the static case proceeded first through 
the introduction of graphical methods by Mayer— and othera^^ in the 1940's and then greatly 
profited from graphical resummation techniques'^ which were subsequently supplemented 
by functional methods'^ as discussed below. Of particular interest here is the functional 
formulation of Percus^"^ which cleanly connects with the widely applied self-consistent 
approximations named Fercus-Yevick-^S and hypernetted chain^ which fit prominently into 
the tool kit of anyone studying the statics of fluids. The theory, which has been applied 
to a large variety of systems, is in the form of nonlinear integral equations connecting the 
radial (pair) distribution function, the direct correlation function and the Ornstein-Zernike 
relation^. 

Why has it taken 40 years to extend the static theory to the dynamic regime? 

Part of the answer is connected with the difficulties in developing an efficient field 
theoretical^ description for classical field dynamics. It was not until the work of Martin, 
Siggia and Rose22, introducing field doubling via conjugate response fields, and its general- 
ization by Janssen^i^ De Dominicis^, Grahan>2^ and others^I, that we had a self-consistent 
treatment for field-theoretic models like the time-dependent Ginzburg Landau models^ and 
all of the models representing dynamic critical phenomena universality classes^. In orga- 
nizing these dynamical theories it is important to carefully incorporate causality. This is 
related to the issue of the proper treatment^- of the Jacobian of the transformation from a 
Langevin description in terms of noise to the path integral description in terms of physical 
fields. There are some apparent ambiguities in determining the Jacobian. This Jacobian 
was identified at least qualitatively in the early work of Onsager and Machlop^. The treat- 
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ment of the Jacobian has led to a fascinating set of extensions of the theory to include 
the topics such as ghost fermions^, supersymmetry^'^ , Onsager's Reciprocity relations^ and 
Jarzynski/fluctuation theorems^. We intend to return to these topics in future work since 
they can be explored in the case of particle models of the type studied here. Kinetic theory 
is complicated^^ compared to conventional field theories because the collective variables, 
particle density and phase-space density, are distributions (sums of 5-functions) not smooth 
fields. This leads to nonlinear constraints like 

p(xi)p(x2) = p(xi)(5(xi — X2) + two particle terms. (1) 

One important aspect of the kinetic theory problem is that static development in terms 
of the density is strongly nongaussian^^. By this we mean that the density fluctuations of 
an ideal gas are not gaussian and one has a set of irreducible vertex functions which can 
not be treated as small. The connected vertex functions for an ideal gas are not small. We 
return to this important point below. 



II. SMOLUCHOWSKI DYNAMICS 



Let us begin by defining the dynamical system of interest. Consider a system of 
particles with configurations specified by the coordinates Ri which satisfy the equations of 
motion 

Ri = DF, + T], (2) 

where the particles experience force 
with total potential 

u{R) = ^Y.yiR^-Rj) (4) 

where we choose ^(0) = V'{0) = V"{0) = and we have suppressed vector labels to 
unclutter the equations. There is a noise source rji for each coordinate which is taken to be 
gaussian with variance 

ivMVjit')) = 2kBTD6{t - t%, (5) 
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where D is a, diffusion coefficient. It is conventional to develop kinetic theory in terms of 
the phase-space density /(I) = J^iLi^i^i — Ri{ti))^{Pi — Pi{ti)) its cumulants. For 
the case of Smoluchowski dynamics this would suggest using the particle density p(l) = 
Y^iLi^i^i ~ Rii'ti)) as the Martin-Siggia-Rose (MSR) field. The key to our development 
is that we break with tradition and treat the particle coordinates as our MSR fields with 
accompanying conjugate response fields. 

In the approach developed here we keep track of degrees of freedom by coupling to them 
with external source fields. In principle we can keep track of all the degrees of freedom in the 
system. In practice, however, we are interested in following a small set of collective variables 
we label $. For this system, the density p is essential since it governs the static equilibrium 
behavior and, from the point of view developed here, is always included in $ = (p, . . .). For 
reasons that will be developed below we must also include in <l> a field B which is constructed 
and interpreted below. The set of collective variables treated = (p, S, . . .)), is fiexible 
and controlled by pairing each observable with a conjugate external field {H — Hp, Hb, ■ ■ ■)■ 
The set $ must include the fields p and B since we need both to control and manipulate 
the interactions in the system. One can include other observables in the problem, like the 
potential energy density, but they play a more passive role in the development. 

We take advantage of the fact that, while the density is strongly nongaussian, the positions 
Ri, may be more profitably thought of as gaussian variables. Therefore, in developing 
our theory we do not work in the Langevin description, but instead work in the MSR 
representation. (As discussed in Appendix A, one has at least three approaches to choose 
from: Langevin, Fokker-Planck and Martin-Siggia-Rose.) The generator of cumulants in the 
MSR representation is given by VFiv[-f^] which is related to the iV-particle partition function 

by 



i=l 



V{Ri)V{Ri)d''Rf 



Po(i?f )e-^«e^-* (6) 



where we have a probability distribution PQ{Rf^) governing the sys- 
tem at the initial time to. The shorthand notation H ■ $ means 

Jd'''Xij^dti[Hf,{xi,ti)p{xi,ti) + HB{xi,ti)B{xi,ti) + ...]. In most of our discussion 
here we assume the system is in equilibrium initially and the initial distribution for a set of 
N particles is canonical: 

Po[i?o] = e-'^^(^°)/Zo (7) 
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where U is the potential energy defined by Eq.(|3]) and (3 is the inverse temperature. The 
MSR action for the problem is given by 

Ar = dtiJ2 \Ri{ti)kBTDRi{ti) + iR,{ti) ■ (Riih) - DFiih))] + Aj (8) 

J to L 

where the contribution to the action Aj is from the notorious Jacobian^. (The steps leading 
from the Langevin description to the MSR field theory description are discussed in Appendix 
lAl ) The Jacobian plays a crucial role in this kinetic problem and is defined by 



6 (dRAf) 



—6(t - t') - \ { 6(t - t') 
dt' ^ ' dRi{t') ^ ' 



(9) 



Exponentiating to write the Jacobian as a contribution to the action givea^ 



where 6^(0) = 1/2. Together, Eqs.([6]), ([8]) and (fTOj) define the problem of interest. 

Now we want to make a separation of the degrees of freedom into two groups; one group 
consists of some collective variables to be chosen, and the second group consists of all the 
rest of the degrees of freedom. The first step in this separation is to rewrite Eq. ( ITOj) in terms 
of the particle density, p(l) = 5{xi — Ri{ti)). We find 

/•oo /• 

Aj = 6(0) J dt J d'^xd'^yDVlpix, t)V{x - y)p{y, t) + constant 

roc r 

= dt d'^xd%B J {x,t)V{x -y)p{y,t) + constant (11) 

Jto J 

where we have defined the quantity 

Bj{x, t) = e{0)DVlp{x, t) (12) 



and the constant can be absorbed into the normalization of the partition function. 
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Next, notice that the dynamic part of the interaction contribution to the action can be 
rewritten in the form 

roo ^ ^ roc r 

D dtY,iR^{t)F,{t) = dt d''xd''yBoix,t)V{x-y)p{y,t), (13) 

where 

N 

Boix, t) = Dj2 ^Ht) ■ ^rA^ - Mt))- (14) 

We can then combine this contribution to the action with the contribution from the Jacobian 
to obtain the dynamic part of the interaction in the form 

Ai = J dld2B{l)V{12)p{2) (15) 

where -B(l) = -Bo(l) + -Bj(l) is the field discussed above, J dl = dtid'^xi, and 

V{12) = V{x,-X2)5{ti-t2). (16) 

Writing things out exphcitly, the conjugate field is given by 

N 

5(1) = DY^ [{RrtVi + ^^(0) V?)] 5{xi - R,{h)). (17) 

i=l 

We can then write the partition function given by Eq.([6]) in the form 

N r 



n 

1=1 



V{Ri)V{Ri)d''R\ 



d d(0) 



Po(i?f^)e-^°-^^+^* (18) 



where Aq is the quadratic part of the action excluding the quadratic contribution to the 
initial probability distribution 

Ao = dhJ2 \RikBTDR, + iR^ ■ R^ . (19) 

Jto L J 

Notice that we have constructed things such that the coordinates are constrained to have 
the values Rf^ at t = to- We then average over these values. Here we are explicitly treating 
the case where the system is in equilibrium at t = to? but more general situations are clearly 
compatible with the development. The interaction part of the action (including the initial 
probability distribution) is given in the compact form 

^i = \T.I dld2^^{l)a^,{l2)^,{2) (20) 
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where the Greek labels range over p and B and we introduce the interaction matrix 

a«,,(12) = (-/3F(12)) [p^p,5{h - to) - {paB, + B^p,)] (21) 
where we have introduced the useful notation 

Pa = (5a,p (22) 

and 

Ba = 5a,B- (23) 
The canonical partition function can be written in the convenient form 

= Tr(^)e-^^+^* (24) 
where we have introduced the average 

Tr(^)C = / n V{R,)V{Ri)(fRf'^ Po{Rf\-^'0{R). (25) 

Notice that the single-particle contribution to the action, Aq, (Eq.( |T9l) ). is included in the 
weight in Tr. Thus the class of problems of interest are defined in terms of a path-integral 
formulation. 

Note that the case where there is a strong external potential acting on the system of 
particles is easily treated within the development. Suppose that the total force acting in 
Eq-® on particle i is of the form 

Fj = F, + (26) 

and the external force is generated by a potential 

Ff = -VrJ d''xUE{x,t)p{x,t) = -VRMMt),t). (27) 

An important practical example is the case of optical tweezers where this external potential 
or trap can be taken to be of the form 

UEix,t) = ^K{t){ix-Roit)f) (28) 
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where K{t) is a controllable amplitude for the potential and -Ro(^) is the position of the trap. 
If one follows the development of the previous section, one finds that the external force 
generates a term in the action 

Au, = J dlB{l)UE{l) (29) 

which can be included in the one-body term in the action. In this case, the initial conditions 
can be infiuenced in several ways. If the trap is turned on for times t > to there is no change 
in initial conditions, while one could prepare the system in a static trap where one would 
need to add a term to the initial potential energy, 

U -^U + J d'^xUEix, to)pix), (30) 

and UE{x,t) = UE{x,to) for t < to. There are many other possibilities. The net result of 
introducing this external potential Ue that couples to the density is to physically produce 
the dynamic coupling to the field -B(l), 

Hb{1) = Ue{1). (31) 

If one is interested in fiuctuations in equilibrium in the presence of a time-independent 
inhomogeneous potential u{xi), then one makes the replacements Hp{l) = u{xi)6{ti — to) 
and Hb{1) = u{xi). Hp adjusts the initial condition and Hb has the effect of changing the 
equation of motion to include the force due to u{x). 

We have succeeded in writing our non-equilibrium problem as a path-integral character- 
ized by a field-dependent partition function written in the very compact symmetrical form 
in the grand canonical ensemble, 

oo N 

Zt[H] = ^ /!^2-^We/'ii^^(i)-*(i)ei/'iW2*(i)-a.$(2)^ ^32) 

which emphasizes the role of the collective fields p and B. We have yet to show that this 
can be expressed in a form which produces a self-consistent form of perturbation theory. 

III. SELF-CONSISTENT DEVELOPMENT FOR THE GENERATING 

FUNCTIONAL 

We now want to rewrite the partition function in a form that allows us to formally carry 
out the average in Eq. fl32|) . We can use the functional identity^ 

^-Aj+H-'^ = e^^e^* (33) 
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where we define the operators 

At = IJ dld2j2^ais{l2)Ha{l)Hi3{2) (34) 



and 

The interaction matrix aai3 is given by Eq. (!2T]) . Then, using Eq. (!33|l in Eg. (1241) gives 

Zj^ = e^^Tr(^)e^-*. (36) 
Next, we restrict the set of fields $j to those that are one-particle additive, 

N 

= E € (37) 

a=l 

which is true for the particle density p and conjugate field B, and notice that the sum 
over the degrees of freedom in Eq. fl36l) factorizes into a product of sums over the degrees of 
freedom of each particle. Together, these observations lead to the result 

ZP = Tr(^)e^* = (Zi)^ (38) 

where the noninteracting partition function for a single particle is 

= rr^e^-^'^ (39) 

Working in the grand canonical ensemble, the grand partition function for the interacting 
problem is given by 

CO y oo , 

- 2. Po ^ - 2. Z^ 

= ^Ar ^ P0_^7V _ ^Ar^Wo (40) 
N=0 ■ 

where po is the fugacity or bare density and 

Wo = poTr(i) e^-<^''' = fr e^'^ (41) 

where in the last line we have dropped the particle label on the trace and the field 0. The 
cumulants of the fields are generated by taking functional derivatives of the generating 
functional 

W[H] = In Zt (42) 
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with respect to Hi. The one-point average in a field is given by 

where we have used a compact notation where i labels space, time and fields p or B. Sub- 
stituting for using Eq. fHU]) we find 

Zt oHi 
Zt 

= 1 rr0ie^^e'^-^e-^^e^^e^°. (43) 
Zt 

It is not difficult to prove the functional identity: 
where 

^=^E^yM- (45) 
is a self-interaction contribution and 

will play an important role as we go along. Using Eq. fHij) back in Eq. (H3|l . we obtain 

Zt 

= f r0,e^-'^+^+^'^[^] (47) 

where 

AW[H] = W[H + F]- W[H]. (48) 



An interesting check on the theory is to show that the self-interaction E given by Eq. ( 14511 
vanishes. This follows if the potential is constructed to be zero at the origin (^(0) = 0). 
Our most important result is given by: 

G, = Tr0,e^-^+^^[^l. (49) 
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Another result useful in this description follows from taking the derivative of Eq. (HO|) 
with respect to po, following steps similar to those leading to Eq.f l49|) . and integrating with 
respect to po which leads to the result 

W[H, Po] = r rfa;Tre^-*+^[^+^'^]-^[^'^]. (50) 



It takes a little calculus to show that the derivative of Eq.f lSOp with respect to H leads to 
Eq.( H9l) . Eq.( l50l) can be rewritten in terms of the more fundamental identity 

■^Zt[H, Po] = Tre^'-^ZTiH + F, po]. (51) 
opo 

What about response functions? The response of the density to an external potential Ue 
which couples to the density is given by 

We must now work to show why Eq.( H9|) is very desirable. There are several ways one 
can use Eq.f l49p to build an approximate theory. In comparison with the static theory one 
would guess that density expansions would be the most successful. This may be so, but 
working with expansions in the pair interaction are conceptionally simpler and more direct. 
It seems clear that in developing density expansions one will be able to make contact with 
the hypernetted chain and Percus-Yevick approximations. There appears much one can do 
about coupling constant renormalization. We return to discuss density expansions elsewhere. 

The dependence of the theory on the pair potential is controlled by the quantity AVr[i7] = 
W[H + F] — W[H]. We can expose the dependence on the potential by constructing the 
functional Taylor-series expansion 

mm = E F.^wm + E \P'F,^nH] + . . . (53) 

and we can conveniently introduce the set of cumulants: 

G^.-^ = ^^---4rmH] (54) 



to obtain 



5Hi 6Hj SHk 



AW[H] = ^ F,Gi + J2 \FiF,G,, + ^ j^F^F^Fj^G^^u + • • • (55) 

i ij ijk 
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with Fi given by Eq. (H6|) . Clearly, in this form we can take AW to be a functional of G,. 
One can then use functional differentiation to express higher order cumulants in terms of the 
one- and two-point correlation functions Gi and Gij. One has, for example, the manipulation 
expressing the three-point cumulant in terms of lower order objects, 

Gijk J- jj Gij ^ ^ Gi^GjnGkp^mnpj (^^) 
"-'^fc mnp 

where the irreducible three-point vertex is given as a functional derivative of the two-point 
irreducible vertex 

and Tij is precisely the matrix inverse of the two-point cumulant 

^^^ikGkj = Sij. (58) 

k 

The beauty of the modern field theoretical development is that an approximation for the 
two-point vertex as a functional of the Gi and Gij generates self-consistent approximations 
for all higher order correlation functions. The method is set up to carry out various types 
of renormalization like replacing the bare interactions with effective interactions. This will 
be exploited elsewhere. We expect the situation here to be similar to quantum many- 
body theory. Self-consistency and conservation laws can be brought together to suggest 
ways of generating approximations as done by Kadanoff and Baym with their <l>-derivable^ 
approximations. A key constraint is equilibrium is the fluctuation-dissipation theorem. 

Some of the structure of the theory can be appreciated via the establishment of a dynamic 
generalization of the static Ornstein-Zernike relation^. Starting with the functional equation 
for the two-point cumulant, one can use the chain-rule for functional differentiation to obtain: 

k \6Gk ) dHj 



k 



(59) 



14 



where 

g,, = frMje"-"^^^"^ (60) 

is a single-particle quantity and we have the memory function^, self-energy, or dynamic 
direct correlation function given by 

= f r0,e^-^+^^-^AIV. (61) 

Since AW can be treated as a functional of Gi we see at this stage that we have available a 
self-consistent theory. If we define the matrix-inverse 

likGkj = Sij (62) 

k 

then the two-point vertex is given without approximation as 

^ij = lij + Kij (63) 

where 

Kij = -^'-fikCkj. (64) 

k 

IV. NONINTERACTING SMOLUCHOWSKI SYSTEM 

The first step in applying this theory is to work out the noninteracting cumulants for 
the fields $ = {p,B). This calculation for the fundamental objects R, R is carried out in 
Appendix [Bl The cumulants for the collective fields are worked out in detail in Appendix 

o 

The final results are 

^ i=l 

where 



GV...Bp...p{^. ...,£,£+ 1, n) = po(27r)^5 9^ m ■ ■ ■ K^)^ " (65) 



and 



n n 

Nn=-DYll^^■^^i^-h\ (66) 
^ i=l j=l 



b{j) = D Y 1i- - t,) (67) 
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and where D = ksTD. Explicitly, the two-point cumulants are 



Gf^{q,q'-tX)=P,{2irY5{q + q')e-^^' 



t-t' 



(68) 



Gfsiq, t') = -poi2n)%q + qWO{t - tOe'^^'^*"*'), 



(69) 



G'^i{q,q';t,t') 



Poi2nf6iq + q')Dq^e{t' - t)e'^^"^''-'^ 



(70) 



and 



{q,q';t,t')=0. 



(71) 



Notice that GpB is retarded and proportional to g^, while Gb...b = 0. The results for 
density cumulants agree with the results from recent work— that shows statistical dynamics 
of the density of noninteracting Brownian particles can be described by a cubic field theory 
where the density is the fundamental field. 

V. PERTURBATION THEORY FOR THE TWO-POINT CUMULANT 

The perturbation theory can be organized in terms of the irreducible vertex functions. 
It is clear from the generalized Ornstein-Zernike equation, Eg. (1591). that the matrix inverse 
of the two-point cumulant is given by Eg. (1631) with the matrix 7 defined by Eg. (1621) and 
the self-energy Kij by Eg.llMl). To get started, one constructs the noninteracting two-point 
cumulant using Eg. fl65|) and one finds the matrix inverses to be given by 




(72) 




(73) 




poDkf 



(74) 



(75) 
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Working to first order in zero external field, AW = J2u PuGu and one has contributions 
to Qij given by 

G,j=fr(t)43{l + ^W + ...) 
= fr0,0,-(l + ^F„G„) 

u 

= G^f + Glff^akuGu- 

k,u 

It takes some manipulation to carry out the various contributions, as will be discussed in 
detail elsewhere, but ultimately 

Po 

where p = (p) is the physical average density, corrected at first order to be p = po/(l + 
poPV{q = 0)). This follows from the perturbation theory analysis of Gj. 
Turning to the dynamic direct correlation function, we have at first order 

= Fj = Y,crjkh- (76) 

k 

Putting this result directly into the defining equation for the dynamic direct correlation 
function, Eq. flGl]) . gives 

c\f = frcpi E (^jk^k = E G^l^akj . (77) 

k k 

The contribution to the two-point irreducible vertex is given by the very simple result 

K^^' = -T.ll^'G^a,, = -a,,. (78) 

k,e 

Let us look at the first-order theory for the two-point correlation function. It satisfies, 
where it is understood that po is replaced by p, the matrix kinetic equation 

Gij = G^f + E Gfi}akiGij. (79) 

kt 
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This is a matrix equation which holds for times ti,tj > to and the two-point cumulant is 
written more exphcitly as 



where a, takes on the values p and B and the translational invariance of the system is 
reflected in the multiplying ^-function. 

Traditionally, there have been two recipes or protocols^'^ for evaluating the two-point 
cumulant. Both reduce the problem to effectively a one-time problem. The kinetic theory 
protocol (KTP) is to treat the problem as an initial value problem with the system in 
equilibrium at time to and we determine the single-time correlation function Gpp{q,ti — to)- 
This quantity is available in the current approach by setting t2 = to in Eq.(l59l) and using 
Laplace transforms. Traditionally, one organizes kinetic theory via the time correlation 
function 



where L is the Liouville operator-^ in the case of Newtonian dynamics. 

In the second protocol, called the field theory protocol (FTP), one takes to — — cxo 
and builds up the equilibrium structure from the noise. One of the technical advantages 
of this approach is that one can maintain time translational invariance over the time line 
and there is only one time in the problem, ti — t2. It is natural to work in terms of time 
Fourier transforms in this case. This allows one to understand the causal structure in terms 
of properties on the complex plane. Our theory here is similar to quantum many body 
theory where one builds up the equilibrium correlation using thermal Green's functions^. A 
difference is that in the quantum case in equilibrium one must satisfy the KMS boundary^ 
conditions. Both protocols are included in the development here. It offers the opportunity 
of developing approximations that are internally self-consistent and one would prefer both 
procedures to produce the same results. 

One can work out the full solution to the set of matrix equations given by Eq.( l79|) with 
the simple result 



(80) 



C^B(g,t) = (5V''*A) 



(81) 



Gpp{q,h,t2) = S{q)F{q,\ti-h\) 



(82) 



where 



F{q,t) 



(83) 
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The static structure factor— is given by 



S{q) 



P 



(84) 



(l + p/3r(g)) 



and the physical wavenumber dependent diffusion coefficient is given by— 



D{q) = Dp-pS~\q). 



(85) 



Notice that we can, at this level of approximation, introduce the notion of an effective 
potential. Comparing Eg. (18^ with the static Ornstein-Zernike^ relation we can identify 
the effective interaction 



where cj:,{q) is the physical direct correlation function which is assumed to be known by 
other means. We can, for example, assume that C£,{q) is given in the Percus-Yevick approx- 
imation for hard spheres^. With this effective interaction one can work out the results of 
perturbation theory in Veff{.(i)- 

One can also use the two protocols discussed above to analyze Eq.f l79p . In the KTP 
one sets t2 = to and notices that only retarded quantities remain in the kinetic equation 
which can be solved directly by taking the Laplace transform. In the FTP where one takes 
to — — oo, ones sees, after taking the limit, all reference to the equilibrium static structure 
is gone and one has time translational invariance. After Fourier transforming over time, the 
equations are reduced to a set of algebraic equations which are simply inverted to give the 
same solution in the frequency regime. Inverting the Fourier transform leads back to results 
found from the complete two-time solution. 

To demonstrate the versatility of the method, consider a system initially (t = to) in 
equilibrium at temperature Tj, but with noise driving the system at T for t > tg. At 
first order in perturbation theory, the system still satisfies Eq.( !79|) but with the t = tg 
contribution in a at temperature T. The solution of this problem is only slightly more 
complicated than the equilibrium case. One finally has the solution 



Gpp{q, ti, t2) = S{q)F{q, |ti - tsl) + {So{q) - S{q))F{q, h - to)F(g, ta - to) (87) 



where So{q) is the static structure factor at temperature T/ and S{q) is the static structure 
factor at temperature T. Now time-translational invariance is broken, but is restored as the 
system decays to equilibrium at temperature T. 



VEFpiq) 



r'coiq) 



(86) 
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It was claimed earlier that this method could provide self-consistent approximations for 
higher-order cumulants. Within the first-order theory one can generate expressions for the 
triplet correlation functions. We easily find for the first order theory that 

vup 

where these are the first order Gs on the right hand side and the zeroth order three-point 
vertex is given by 

Ivup ^ y liv 1 ju Ikp ^ vup C^'^) 

vup 

where all zeroth order quantities can be evaluated using Eq.( l65|) and Eqs.( !72|) - (175|) . 

At second order in the effective potential we have two contributions to the two-point 
vertex. The first piece comes from the self-energy contribution to the dynamic Ornstein- 
Zernike equation, Eq.( l6T]) . where, keeping the second order terms, we have 




— — X! ^i^i X! C^ii^ikpGpj (90) 

ij e,p 

and, in the simplest second order approximation, we replace the 3-point vertex with the 
zeroth order result. After some simple manipulations that will be described in detail else- 
where, the second order contributions to the two-point vertex coming from the self-energy 
can be written in the mode coupling form 

= E llklSGkp6G,4Z (91) 

kpin 

where the new element is that the high wavenumber convergence of the integrals comes not 
from the three point vertices but from the subtraction inH 

SGkp = Gkp — G^j^p . (92) 

In some field theoretical treatments'^ of the DK model one obtains memory function kernels 
which do not vanish in the noninteracting limit. To complete the second-order model we 
must also work out the second-order terms coming from Q. This term is linear in the full 
two-point correlation function and if there is an ENE transition at this order, it could show 
the stretching phenomena found phenomenologically in Goetze's F12 model'^. The complete 
second order model will be treated in a separate pubhcation. 
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VI. CONCLUSIONS 

We have presented here a reformulation of kinetic theory which is self-consistent. It allows 
one to study problems which are difficult to treat using other methods. 

1. We outlined a clean derivation of the mode coupling model at second order in pertur- 
bation theory. We will analyze whether this model supports ENE transitions elsewhere. 

2. The method presented here allows for a systematic method for analyzing corrections 
to this second order result including higher order correlation functions. 

3. The method allows one to treat nonequilibrium problems like temperature quenches 
as demonstrated above. 

4. Completely unexplored is the fact that the perturbation theory has been developed in 
the presence of space and time dependent external fields. Thus this method could be useful 
in problems of optical pinning and highly inhomogeneous situations. 

5. It seems likely that this method will be useful in treating meta- and unstable systems. 
A first approach within perturbation theory is to formulate a dynamical van der Waals 
theory. Similarly it seems likely that this approach will be useful in developing a dynamic 
theory of melting once a few ideas from density-functional theory are integrated into the 
development. 

6. Our focus here, because of its simplicity, has been on Smoluchowski dynamics, but as 
will be discussed elsewhere, the method developed here can be applied both to Newtonian 
dynamics and Fokker-Planck dynamics as well as a broader class of models^. In the case of 
Newtonian dynamics, this approach offers an alternative to the conventional development in 
terms of the Liouville operator. It may be important that in this approach it is not necessary 
to use basis states in perturbation theory labeled by a continuous momentum index. 

7. It is clear that we can study using these methods the mapping of Fokker-Planck dy- 
namics onto Smoluchowski dynamics in the large mass limit. This is of course an elaborated 
version of the calculation leading to Einstein's relation^ between the friction coefficient 
and the diffusion coefficient. Of particular interest to us is whether in integrating out of 
the momentum degrees of freedom one generates density dependent diffusion coefficients. 
The effects of such a density dependence been shown to be physically important within the 
random diffusion model.—. 

8. All of this development is compatible with the bizarre developments initiated with the 
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introduction of ghost fermions^ into the treatment of stochastic dynamics. We expect the 
subsequent developments hke supersymmetry^ to manifest themselves as in field theories 
where one can make a connection with Onsager's reciprocal relations^ and the fluctuation 
theoremsA^ found in the strongly nonequilibrium regime. 

9. Independent of the relevance of point 8, the connection to the fluctuation theorems 
and Jarzynski equalities should be explored carefully. 

10. These methods allow one to study one-point quantities in more complicated situa- 
tions. In the calculation outlined above for homogeneous systems we did not mention the 
self-consistent determination of the equation of state because the results are rather dull. 
This will not be the case for inhomogeneous systems. 

11. It seems clear that a more quantitative method can be established if one determines 
the best way of expanding AW in a density expansion. This would make close contact to the 
work of Percus and the standard approximations of Percus-Yevick and hypernetted chain. 

12. This method allows us to investigate the claim by Das and Mazenko^ that momenta 
are associated with a mechanism which cuts off any possible ENE transition in conventional 
fluids. 

13. Even in the case of Smoluchowski dynamics, the connection between the theory at 
the microscopic level and fluctuating nonlinear hydrodynamics has not been established. 



Can the noninteracting field theory of Ref. 



52 



be extended to the interacting regime and 



connected up to the microscopic theory treated here? 

14. The existence of fluctuation-dissipation theorems is very important in equilibrium. 
This will be explored in depth elsewhere. 
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Appendix A: Langevin, Fokker-Planck and MSR Dynamics 



1. Langevin Model 



The Langevin model we study is given by 

dRi{t) 



dt 



DFi{t)+r]i{t) 



which is vahd for t > to and has the initial condition 

Ri{to) = 



(Al) 



(A2) 



governed by the probability distribution Po[Rf^]. Fi{t) is the force acting on Ri{t) and is a 
local function of Ri{t). We assume that the noise T]i{t) is gaussian white noise with variance 



{r],it)r],{t')) = 2kBTD5i,5it~t') 
where D is a diffusion coefficient. The associated probability distribution is given by 



P[rj] = Noexp 



dt 



(A3) 



(A4) 



Ho AD 

where Nq is a normalization factor. The partition function in an external field is given by 



Z[h] = I Vr]P[7]] J I?i?(°)Po[i?^°^]5(/i ■ R) 



(A5) 



where the time-dependent external field couples to the system via 



S{h ■ R) = exp 



N 



dtJ2Ri{t)hi{t) 



(A6) 



and where, via the Langevin equation and its initial condition, the field Ri is a functional of 
the noise and the initial condition Rf*^ . Our goal is to determine the generator of cumulants, 
W[h] = In Z[h]. 



2. Fokker-Planck Dynamics 

It is convenient to analyze this set of dynamical models using the Fokker-Planck descrip- 
tion. We follow here the development of Kim and Mazenko^-. If we define the field 

N 

gdt) = I[H<P^-R^it)), (A7) 
i=l 
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we may then consider the time correlation functions 



(A8) 



where the average is over the noise and the initial condition. 

Taking the time derivative of G^^^'{t) using the chain-rule for differentiation and the 
Langevin equation, Eg. flAip . one is left with 



It is then not difficult to show, remembering that rj is gaussian, that 



(A9) 



(AlO) 



where we have used D = ksTD, the result 



SRiit) 1^ 

— -Oi 



(All) 



Hit) 2^ 

which follows from Eg. flAip and the assumption that the initial field configuration is inde- 
pendent of the noise for t > 0, 

SRM 



Hit) 



0. 



Using Eg. flAlOp back in Eg. flA9p we can write 

d 



dt 



G^^(f)'it) — D^G(j)^^'{t), 



where the Fokker-Planck operator is defined by 

6 



The formal solution to Eg. flAlSp is 



6 (pi 



(A12) 



(A13) 



(A14) 



= e^^(*-*°)[5(0-0')^o(0')]- 
Integrating over all cj)' gives the eguilibrium probability distribution 



(A15) 
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e^^(*-*«)Po(0)- (A16) 



It is easy to see that the equihbrium solution is given by 



where 



50, 

Eq. (IA15|) then takes the form 



P^(0) ^W^ = -—- (A17) 



= -1^- (A18) 



G^,^it) = e''^'[Si^-^')W^]. (A19) 
Any two-time correlation can then be written in the form 

CAB{t) = {A{R{t))B{RiO))) 

Vct)V(t>'A{(t>')B{(f)G^^^\t) 

V(l)B{(f))e^^'A{(f))W^ (A20) 

where in the final step we assume the system is in equilibrium. 

In our development it is useful to introduce the adjoint Fokker-Planck operator 



Fi - ksT— 



(A21) 



i 

If the equilibrium average is defined as 

{A) = J V{(j,)W^A{(j>), (A22) 

one can show that 

CABit) = J V{(j))Bi(j))e''^\A{(j))W^) = {B{(j))e^^'Ai(j))). (A23) 
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3. Multiple time correlations 

Consider the multiple time correlation function 

G^0O,01,02v,0n(^O) ^1) ^2; • • • ) ^n) = {9<j> n 

where g^pit) is defined by Eg. flA7p and we assume t„ > > ■ ■ ■ > t2 > ti > Iq. Using the 
same approach as used for treating the two-point quantity we have 



X 



-DFjG0Q_0-^^02^..._<^^(to, ti, t2,..., tn) 



-J2^BT{r]i{tn)g^^_,{tn-i)---g^2it2)g4>i{ti)g^oito))]- (A25) 

i 

Because of causality, g^X^d f*^^ < is independent of the noise at t„: 

{vii^g^iQ ■ ■ •fl'9i2(^2)fl'0i(ti)fl'0o(^o)) 
= D{^£^^g^^_,{tn-i) ■ ■■gUh)g^Ati))gM)- (A26) 

The treatment of is the same as the two-time case and we obtain the result 

d 

^G0o,0-^,02,.-,</'n (^0,^1,^2, ■ ■ ■■,tn) = -D<^„G'<^o,<^j,02,...,0^(ti,t2, ■ ■ ■ ,'tn) (A27) 
dtn 

where is the Fokker-Planck operator given by Eg. flA14l) . This has the formal solution 

G^(^O,</>l,</>2v,0n-l,</'n(^O) ^1) ^2? • • • tn—1, tn) 

= e^^"(*"-*"-^)G'^oA,02,...,</'„-i,</'n(^o, ^1, ^2, • • • tn-U tn-l)- (A28) 

However, 

G</,o,<^l,<^2,...,</,„_i,0„(to,^l,^2, • • -tn-l, tn-l) = {g<t,„{tn-l)g4>„_i(tn-l) ■ ■ " 5'</'2 (^2)fi'</.i (^l)5'</,o (^o)) 

= (^(0n — 0n-l)G<^o,0l,02,-,</'n-i(^O5 ^1) t2, . . . tn-l)- (A29) 
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Clearly we can work this out recursively: 



X - - 0„_2) . . . 5(01 - (f)o)W^g. 



(A30) 



If we introduce the notation 



(A31) 



we can write 



G<|,o,<|,l,<|,2,...,<f>r^-l,<f>Ato,tl,t2, . . .tn-l,tn) — U^^-^^_^{tn - tn-l)t/0„_i;<^„_2 (^n-l " tn-2) 



(A32) 



This is the result we need in developing the path integral approach. 

4. Path Integral Form 

How is the partition function in a field related to these multiple-time correlations? In 
Z[h], given by Eq. OASp . we make the special choice for the external field 



which amounts to dividing up the time interval into a grid. Next, multiply by 



(A33) 



s=0 



1 "l 



and integrate over A*. Then, we have 



•' s=0 ^-8=0 
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s=0 



X 



= / Pr/PM / I)p(°)Po[i?^°^]n^?0.(ts)e^-''^'''^*^^ (A34) 

s=0 

where again the are 5-functions which allow us to make the replacement 

eE,, ''f^^C*^) = e^^>» '^^'^'^ = S[h ■ <j)] (A35) 
which comes out from the average over noise and initial conditions. We have 



'-8=0 



= S[h ■ (f)]G^,,^,,,„^^„{to,ti, . . . ,tn). (A36) 

We finally obtain the result for the partition function we want by doing the functional 
integral over 0: 

Z[h]= I P0^[/i-0]G^„,<^,,..„<^„(to,ti,...,tn). (A37) 

The key here is to notice that we have an explicit expression for the multi-time correlation 
function G</,o,(/)i,...(/)„(^05 ^i, • • • ytn)- Inserting this result into Eq. (IA37p . we obtain 



Z[h] = J S[h ■ 0]f/0„;^„_i(tn - tn-l)f^</,„_i;^„_2(in-l " tn-2) " " " 

X U^rMt2 - ti)U^,,^,(ti - h)W^,y (A38) 

There are no constraints on the choice of time slices. Let us take the slices to be uniformly 
divided, tg+i = + and we work in the limit of small A and large n. This defines the 
continuum limit for the theory. As a check on the development, notice that the normalization 
as /i — )■ is preserved: 

Z[0] = j V(t)U^^.^^_^{tn - tn-l)U^^_y^^^_^{tn-l - tn-2) " " " 

X U^r,4>iit2 - ti)U^,.^,{ti - to)W^, = 1. (A39) 
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The result given by Eq. (lA38P is naturally interpreted as a path-integral. Let us focus on 
the intermediate time quantities 



(A40) 



Using the integral representation for the 5-function we can diagonalize the the Fokker-Planck 
operator: 



^(/'n;0n-i(^n - tn-l) - e 



Then, we can evaluate 
where 



and 



n-1 



,e " e 



(A41) 



(A42) 



(A43) 



where 



A„ = -D(f)l + ic, 



Putting Eq. fK45l) back in Eq. flMSl) . one has 



Z\h\ = I P0P0I)0on(0o)e^'^'*^^ 



where Aji is the standard MSR action in the presence of an external field, 



A 



R 



D(j)\t) + t(f){t) (f){t) - F{(f){t)) +h{t)(j){t)- 



(A44) 



(A45) 



(A46) 



(A47) 
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Appendix B: Gaussian Single-Particle Problem 

The noninteracting correlations for a system driven by Smoluchowski dynamics is gov- 
erned by the MSR action 



Ao = / dt 

Jto 



oo 



R{t)DR{t) + iR{t)R{t) - h{t)R{t) - h{t)R{t) 



(Bl) 



where h{t) and h{t) are the external source fields and D = ksTD. We then have the 
identities that hold in the range < t < oo, 

S_ 
'5R{t)' 

and 

S_ 

'6R{t) 

Let us begin with the initial condition 



I V{R)V{R)d''RoPo[Ro]-^e-^'^ = (B2) 



r 

I V{R)V{R)d''RoP^[Ro]-^e-^' = 0. (B3) 



Po[i?o] =5(i?o-^o). (B4) 



Evaluating the derivatives of ^40, we obtain 

6 

~ dt 



2DG{t)+i^^^h{t). (B5) 



and 



-i^-Kt) (B6) 



_,dG{t) 
where 

G{t) = {R{t))xo (B7) 

and 

G{t) = {Rit))xo (B8) 

where the averages over R{t) and R{t) are in the range to < t. We must now solve these 
equations to obtain the generating functional. 

Using the initial condition R{to) — 0, we find that 

^ /"OO _ POO 

G{t)^-i dth{t)= dTg{T,t)h{T) (B9) 

Jt Jto 
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and 



G{t) = Xo+ f dt[2iDG{t) - ih{t)] 

Jto 

= Xq-i / dih{i) + 2D dt /_ dt'h{t') 

J to J to J t 

roc ^ POO 

Xo+ drg{t,T)h{T) + dTC{t,T)h{T) (BIO) 

Jtn Jtn 



where 



g(t,t') = -ie{t-t') (Bll) 



and 



C{t,t')=2D dt dt'5{t'-t'). (B12) 

Jto Jt 

The generating functional satisfies 

G{t) = (B13) 
^ ' 6h{t) ^ ' 

and 

_ ^lnZo(/i,/i,Xo) 

^^'^ - — m) — ■ ^ ^ 

The generating functional solution to this set of equations - Eqs. flB9p . fIBlOp . fIBlSI) and 
(]B14p - is given by 



lnZo(/i, /i,Xo) = ^JdtJ dt'h{t)C{t,t')h{t') 

+ J dt J dt'h(t)g{t,t')h{t') + j dth{t)XQ. (B15) 
The full generator requires averaging over the initial conditions: 

Z[h,h] = J d'^RQPQ[RQ]e^^-^-^+^-^'^+^-'3-R° 

= e^^-^-^+^ <^ '^ j d'^RoPo[Ro]e''-'^-^. (B16) 
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All of the equilibrium cumulants can be constructed from Eq. (lB12p as 



C{t, t') = 2D di dt'5{t' - t') 

Jto Jt 



2D f dteit' -f) 

Jto 



2De{t - t') f dt + 2De{t' - t) f dt 

Jto Jto 



2De(t - t'){t' - to) + 2De{t' -t){t- to 



(B17) 
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Appendix C: Collective ^-Correlations 

We need to evaluate ^-correlations in the non-interacting case. We start with 

Wo[H] = fre^-^ (CI) 

which generates all (j) = {p, B) correlations. We proceed by reintroducing the microscopic 
sources h and h and treating 

Zo[H, h, h] = fre" 'f'e''-^+^-^ (C2) 

where 

h-R^ dth{t)R{t). (C3) 
Next, we express the <j) in terms of R{t) and R{t) as 

(t)p{l) = e-'^^i^^*!) (C4) 

and 

0b(1) = -D[{h ■ R{1) + e{Q)kl)ct>p{l)]. (C5) 
Let us introduce the operators 

0^(1) = e"'^'^) (C6) 

and 

0b(1) = &(1)0,(1) (C7) 

where 

6(1) = -D 
so we may write 

and 

Zo[H, h, h] = e^-^Zo[h, h] (CIO) 



dh{ti) 



(C8) 
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where 2'o[/i, h] was determined in Appendix [Bl 

Taking functional derivatives, we can determine all of the noninteracting cumulants of 
the complete set of densities (j). We have 

5 5 5 5 

G'ij...ijp...p(l,...^,^ + l,...,n) 



X 



5Hb{1) ■ ■ ■ 5Hb{^) 5H,{i + !)••• 5Hp{n) 



= . . . + 1) . . . 0,(n)e^^-^-V-5-^ j f/'^/?oPo[i?o]e"-^^-^°|^=,=/,=o- (dl) 

All n-point cumulants have n factors of (f)p and I factors of h corresponding to the number 
of B insertions: 

GB...Bp...p{l . . . ^ + 1, . . . , n) = 6(1) . . . 6(£)0p(l) . . . 0p(n)e^"-^-^^-5-^ 



Because the 4>p{j) are translation operators, it is not difficult to show that 

0p(l) . . . ^p{n)F[h{j)] = F[hij) + (C13) 



where 



Ln{j) = ~tJ2ks5{tj-Q. (C14) 



Thus, we have 



^(h+Ln)-ig-Ro I 



\h=h=0 



G'iJ...Bp...p(l,...^,^+l,...,n) 

6(1) . . .6(£)e5('^+^")-^-('^+^")e('^+^")-5 '^y ci'^i?o^o[^o]e^ 

= 6(1) . . . 6(£)e5^"-^-^"e^"-^ '^y" ci'^i?o^o[i?o]e-^"'^^'^1 

= 6(1) . . . 6(£)e^"e^" ■^'-'^y rf'^i?oPo[^o]e^"-^^-^°|;,=o (^15) 



\h=0 
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where we have defined 



2 

Using the definition of b, we similarly find that 



N^ = l-L^-C- Ln. (C16) 



where 



= -D[k,Ln{t)g{i - tj) - e{{])k]] 

n 

= Dkj ■ J2 ksO{ts - tj) - De{<d)k] 



s=l 



= Dk,- hO{ts-t^). (C18) 
In the average over initial conditions, we need 

n 

Ln-tg = -tJ2ks (C19) 

s=l 

and 



d'RoPo[Ro]e''"-''-''' = d'RoPo[Ro]e K^^' 'J = {27Tf6(j2kA (C20) 

s=l 

which enforces translational invariance in space. Finally, putting this all together we have 

G'B...sp...p(l, 1, n) = poKl) • ■ • me""" i2nf6 ( E A;,) . (C21) 

The argument of the exponential contribution can be put into a more symmetric form. 
Starting with iV„ given by Eq. (lC16p and inserting L„ from Eq. (IC14p . we have 

-inn 

Nn = --T.J:hk,CiU,t,). (C22) 
^ i=i j=i 



The zeroth-order correlation function for the Brownian coordinates is given by Eq. (IB17P 

as 

C{t, t') = 2D [e{t - f) {f -to) + e{t' -t){t- to)] (C23) 
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which, at equal times, reduces to 

C{t,t) = 2D{t-to). (C24) 

In the zeroth order density correlation functions we find the quantity 

1 " 

Nn = ~^Y.k^kJC{t„tJ) (C25) 

with the constraint that J^i^i = 0. This quantity should be time translationally invariant. 
To see this, let us first define 

= 2D (^ti + tj - 2[9{ti - tj)tj + 9{tj - ti)ti]^ (C26) 
= 2D(ti[9{ti - tj) - 9{tj - ti)] + tj[9itj - ti) - 9{t, - tj)] 



2D(ti - tj)sgn(ti - tj) 



= 2D\ti-tj\. (C27) 

Notice that this result holds if t — )■ t — to in Eq. flC26p . 

We then have for the argument of the exponential for the p-B correlation functions 

1 " 

Nn = — — ^ kikjC{ti, tj) 



1 ' 2 

ki Cii + kikjCi. 

i=l i^j 



2 



= — - l^kidi + kikj-[Cii + Cjj — Dij] 

\i=l ijtj 

1 / " 2 "2 1 

= — - {^^kiCii — ^^kiCii — kikjDi 

= D]^Y.kik,\ti-tj\. (C28) 

For the special case of n = 2, we have 

N2 = -Bk\\t^-t2\. (C29) 
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By kinetic theory we mean the statistical dynamics of N classical particles. See G. Mazenko, 
Nonequilibrium StatisticalMechanics, Wiley, (2006) for a discussion of kinetic theory in- 
cluding the modern correlation function approach. 

Brownian motion is interpreted here as a dynamics that is noise driven. For an individual 
Brownian particle located at position R acted upon by noise r), one has R = r). This is clearly 
the continuum version of the random walk problem. 

Field theory methods enter via the techniques for organizing sclf-consistcnt perturbation theory 
in terms of the roles of external fields, cumulants and irreducible vertex functions. The fun- 
damental reference is J. Zinn-Justin, Quantum Field Theory and Critical Phenomena, 
fourth edition, Clarendon Press, Oxford (2002). We henceforth refer to this reference as JZJ. 
Kinetic theory's early history is briefly reviewed in G. Mazenko, Equilibrium Statistical 
Mechanics Wiley, (2000). See page 101. 

By self consistency we mean that the correlation and response functions of interest satisfy kinetic 
equations where the collision kernels, self-energies or memory functions can be expressed in 
terms of the full physical correlation and response functions. In the areas of critical dynamics 
and the liquid-glass transition, this is essential in generating symmetry breaking solutions. This 
can generate a multiplicity of solutions not available in bare perturbation theory. 
Boltzmann's Stosszahlansatz is a prime example of a decoupling approximation. In the colli- 
sion integral of the Boltzmann equation he assumed that the two-particle distribution function 
factorizes into a product of two one-particle distributions, /2 — )■ /i/i. 
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weak supercooled states." 
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it conforms to the details of mode coupling theory. The defining property of an ENE tran- 
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(2001) where they comment: 

"Despite its remarkable practical success, the presence of apparently uncontrolled 
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Smoluchowski dynamics has become identified with the over damped kinetics in colloidal systems 
where the momenta become equilibrated much faster than the positions and one has a dynamics 
which is subsequently organized in terms of the positions. This is the point of view taken by 
A. Einstein, Ann. d. Physik 17. 549 (1905) in his seminal treatment of Brownian motion. M. 
V. Smoluchowski, Phys. Zeit. 17, 557 (1916), still working in coordinate space, generalized 
Einstein's work to include external forces acting on the Brownian particles. These results were 
valid for long times. P. Langevin, Comptes. rendus 146, 530 (1908) shortly thereafter developed 
the Langevin equation description for the momenta of the Brownian particle. Working in terms 
of probability distributions in velocity space, A. D. Fokker, Ann. d. Physik 43, 812 (1914) and 
M. Planck, Sitz den Preuss. Akad. 23, 324 (1917), established the simplest form of the Fokker- 
Planck equation. These works were tied together by a treatment of position and momentum by 
O. Klein, Ark. Mat. Astron. Fyr. 16, 5 (1922). This work was repeated by H. Kramers, Physica 
7, 280 (1940). 
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for more than one Brownian particle was developed by J. Deutsch and I. Oppenheim, JCP 54, 
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Important papers treating the many-particle Smoluchowski dynamics system include B. J. Ack- 

erson, JCP 64, 242 (1976); W. Dietrich and I. Peschel, Physica A95, 208 (1979). 

In the case of Smoluchowski dynamics it has been shown by B. Cichocki and W. Hess, Physica 

A141, 475 (1987) that the memory function can be rewritten in terms of an irreducible memory 

function. This argument was generalized by K. Kawasaki, Physica A215. 61 (1995). For our 

purposes here the point is that there is not a unique memory function form at lowest nontrivial 

order in perturbation theory. Thus for a variety of models one can rearrange perturbation theory 

to have model which supports an ENE transition. One must go to higher order to check the 

self-consistency of an ENE transition. 

The theory we develop here has great potential for treating higher order correlation functions. 
The reason why this is interesting is because there has been a significant amount of work 
associated with the concept of dynamic heterogeneity, e.g. S. C. Glotzer, J. Non-Cryst. Solids 
274, 342-355 (2000). Out of this research has come the idea that the order parameter for the 
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to a diverging length scale £ at the ideal glass transition. They find an upper critical dimension of 
six in their calculation. This work is very provocative since at the level of the dynamic structure 
factor it is well known that MCT does not contain a large length as one approaches the ergodic- 
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look at the four-point quantity C4, to find this diverging length. However this calculation paints 
a picture rather than gives the results of a rigorous calculation. One of the goals for the theory 
developed here is to connect the collective behaviors found at the two-point level to the behavior 
of the three- and four-point density correlation functions. 
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Once one has the ghost fields as described in Ref. 
field, one can search for symmetries of the action reflected in the whole space of the enlarged 
system. For a large class of stochastic systems one finds a supersymmetry which mixes fermionic 
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and conventional field degrees of freedom. The key references are given in Ref. 
further in the development as indicated by the title "Ward Takahashi Identities and Fluctuation- 
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